X(3872) and the bound state problem of D°D*° (D°D*°) in a chiral quark model 
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The bound state problem of D°D*° (D°D*°) is relevant to the molecular interpretation of the 
X(3872). We investigated this problem in a chiral quark model by solving the resonating group 
method equation. We found the system is unbound through S-wave tt and a interactions. The 
inclusion of p and uj meson exchanges is helpful to the formation of a molecule. Because the binding 
energy relies on the coupling constants, we cannot draw a definite conclusion whether a molecular 
state exists in D°D*° (D°D*°) system. When moving on to the bottom counterpart, we obtained 
an S-wave BB* state. 

PACS numbers: 12.39.-x, 12.40. Yx, 13.75.Lb 



I. INTRODUCTION 

In recent years, a number of new charmonium-like 
states have been observed [2, 0, S 0| ■ O ne °f the most 
interesting states is the X(3872). The Belle Collabora- 
tion first discovered this state in the ir + Tr~J/ip channel 
of B decay in 2003 @. Thereafter, CDF @, DO 0, and 
BaBar [8J collaborations have confirmed its existence. 
The X(3872) is almost on the threshold of D°D*°. Its 
width is very narrow (T < 2.3 MeV from the Particle 
Data Group @). The measurements from Belle and 
CDF 11] favor the quantum numbers J PC = 1 ++ , but 
2 h have not been ruled out yet. In the search for a 
charged X state, BaBar excluded the isovector hypothe- 
sis |12j |. 

Experiments have accumulated much information 
about the decay of the X(3872). The analysis from 
CDF 13] supports that the two pions in the channel 
-Y(3872) — ► J/ip come from the p meson. In ad- 

dition, Belle observed the 37r decay w + tt~ w° J / 4> and the 
radiative decay yJ/ifi [llj]. BaBar also reported the ev- 
idence of the latter mode The measured ratios in- 
clude Q 



£(X(3872) -> tt+tt-tt J/VO 



B(X(3872) 



■J/1>) 



and 



6[X(3872) -> jj/tp] 
B[X(3872) J/^tt+tt- 



S[Y(3872) -> jj/tp] 
B[X(3872) -> J/V'tt+tt- 



1.0 ±0.4 ±0.3, 



0.14 ±0.05 



0.3. 



(1) 



(2) 



(3) 



One notes the ratio between the 37r mode and the dipion 
mode in Eq. |T]) indicates the large isospin violation when 
the X(3872) decays. 



Recently, Belle announced a new near-threshold en- 
hancement with M = 3875.4 ± 0.7±H MeV in the chan- 
nel B -> X(3875)if -> D D°tt K This state has 
been confirmed by BaBar [l^]. It is unclear whether or 
not these two X states are the same one. If the X(3875) 
is identical to the X(3872), there are two more ratios [17| 



B[X -» D°D°tt } 



B[X 



B[B a 



-JH>] 
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j+3.1 
'-3.6) 
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1.6. 
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The exotic properties of the X(3872) have triggered 
heated discussions about its nature. The low mass, the 
extremely narrow width and the large isospin violation 
decay are difficult to understand in the conventional cc 
assignment [HI . Up to now, there exist many interpreta- 
tions: a molecular state |20|. l2ll . |22[ l23l. l24j . a cusp l25l| an 
S-wave threshold effec t 1261 ] . a hybrid charmonium [27| . 
a four quark state [1^, |20(] . a vector glueball mixed with 
some charmonium components [3fj| and a dynamically 
generated resonance [3l|. In addition, there are discus- 
sions that the puzzles for the X(3872) may possibly be 
resolved in the scheme of mixing 

The molecular interpretation is the most popular one 
in understanding the structure of the X(3872). In fact, 
the existence of a loosely bound molecule (deuson) of two 
heavy mesons has been proposed long ago [13, [H, 
In such systems, the contribution from the kinetic term 
is lowered because of the presence of the heavy quarks. 
Since the attraction from the light quarks is unaffected 
by the mass of the heavy quark, the formation of the 
heavy deuson is possible. According to the calculation in 
Ref . [3!| , several deusons of two bottom mesons should 
exist while the predicted deusons of two charmed mesons 
are close to the thresholds. 

The proximity of the mass of the X(3872) to the thresh- 
old of D°D*° (D°D*°) motivated its molecular interpre- 
tation. Numerous discussions have tak en p lace within 
this picture 0, EH H S, El E H S3, Si- The mass, 
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the quantum numbers J PC , the isospin violating decay 
and the 3-7T decay appear to be naturally understood. 

However, the ratios in Eq. ^ and Eq. © challenged 
the molecular interpretation. Both are inconsistent with 
the molecular picture's prediction which is around 7 x 
1(T 3 . If the X(3875) and the X(3872) are the same state, 
the values in Eq. (j4|) and Eq. |5j) are also much larger 



than the theoretical predictions. 



B[x- 

B[X- 
B\B° 



T+7T .// 1p\ 

^XK 



TXK+\ 



from the 
is less than 



molecular assumption is 0.05 and gfgr 

0.1 [H. 

Therefore, whether the molecular picture is correct or 
not remains inconclusive. This question is relevant to 
whether D°D*° (D°D*°) can form a molecule. Up till 
now, the dynamical studies of this system are still scarce. 
In the calculation of Swanson [23[ and Wong [22j, bind- 
ing is possible when the short-range quark-gluon force is 
considered. However, the purely molecular assumption 
of the X(3872) was questioned in Ref. [33j |. 

In order to further understand the nature of the 
X(3872), it is worthwhile to study dynamically the molec- 
ular assumption for the X(3872) with various methods. 
In a previous work [4^] , we have investigated at hadronic 
level whether the formation of a bound state of D° and 
D*° is possible. We found that one pion and one sigma 
exchange interactions could not bind the system to an 
S-wave molecule. The same framework was also applied 
to the newly observed Z+(4430) [5(J. 

In this paper, we reanalyze whether the X(3872) could 
be an S-wave D°D*° molecule in a different approach. 
We will study this system in a chiral constituent quark 
model and calculate the binding energy by solving the 
resonating group method (RGM) equation [fO], [52j ■ 

The chiral quark model [53j is a useful tool in con- 
necting QCD theory and the experimental observables, 
especially for the light quark systems. This phenomeno- 
logical model has been quite successful in reproducing 
the energies of the baryon ground states, the binding en- 
ergy of the deuteron, the nucleon-nucleon (NN) scat- 
tering phases and the hyperon-nucleon (YN) cross sec- 
tions. In this model, the interacting potentials between 
the two constituent quarks include the confinement, the 
one-gluon exchange (OGE) part and the pseudoscalar 
and scalar meson exchange part. It has been controver- 
sial whether OGE or vector-meson exchange dominates 
the short-range quark-quark interaction in the low-lying 
baryon states. Thus the vector meson exchange part has 
been included in Ref. [54J. The model was named as 
the extended chiral SU(3) quark model. It was found 
that the OGE is nearly replaced by the vector meson ex- 
changes. By solving the RGM equation, the experimental 
observables were well reproduced. 

Recently, the chiral quark model has been extended to 
study bound state problems for the baryon-meson sys- 
tem (55| and baryon-antibaryon system [56f by solving 
the RGM equation. In this work, we will study similar 
problem for the D°D*° system within this approach. 

This paper is organized as follows. In Section [TTJ we 



present the formalism for the calculation. In Section HTT1 
we give the methods to determine the parameters and 
their values. Then in Section IIVI we show the numer- 
ical results, and the last section gives a summary and 
discussion. 



II. FORMALISM 

A. The molecular picture 

The heavy molecular state bound by the one-meson 
exchange interaction in the chiral quark model can be 
depicted in Fig[T] where A and B are two heavy mesons. 
The OGE and the confinement interactions occur inside 
the color-singlet mesons only. The interactions between 
the two clusters are induced by the one-meson exchange 
potential between light quarks. 




FIG. 1: The molecule bound by one meson exchange interac- 
tion in the chiral quark model. The solid (empty) dot is the 
heavy (light) quark or antiquark. 

If the X(3872) is really a molecule, the wave function 
in flavor space should be @, l49l | 



X(3872) = 



«o 



V2 
y/2 



D°D*° ~ D*°D° 



D + D *- _ £)*+£>- 



(0) 



where the ellipsis denote other hadronic components. Be- 
cause of the large isospin violation in the decay, one ex- 
pects the first part dominates with do >> a\. In the 
following calculation, we study whether D° and D*° may 
form an S-wave molecule with the flavor wave function 

m 



\Xd) = 



1 r 



V2 



\D°D*°) - \D* a D°) 



(7) 



If the answer is yes, this molecular state should lie 
below the threshold and identifying the X(3872) as an 
X^-dominated molecule is favored. Otherwise, the pure 
molecular interpretation of the X(3872) is problematic. 
We search for an answer by calculating the binding en- 
ergy of the system D°D*°. Besides the pion and sigma 
exchange interactions, the p and the uj exchange effects 
are also considered and discussed. 
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Hamiltonian 



The details of the chiral SU(3) quark model can be 
found in Refs. [H, HU. Here we just present essential 
constituents for the calculation. The Hamiltonian has 
the form 



H = 



i=i 



T l -T G + V OGE + V conf + yM ( 8 ) 



M 



where Tj is the kinetic term of the zth quark or antiquark 
and Tg is the kinetic energy operator of the center of 
mass motion. 

The potential of the OGE part reads 



Where Gm is the G-parity of the exchanged meson and 
9ch A2 ™ 



C(g ch ,m,A) 



4tt A 2 — m 2 
A. 



Xi(m,A,r) = Y(mr) Y(Ar), 



(15) 
(16) 



X 2 (m,A,r) = Y(mr)-[ — ) Y(Ar), (17) 

m , 



Y(x) = — . 



(18) 



The tensor term and the spin-orbital term have been 
omitted in the potentials since we consider only S-wave 
interactions. We use the same cutoff A for various 
mesons. Its value is around the scale of chiral symme- 
try breaking (~1 GeV). 



V, 



OGE 

qQ 



g q g Q F q • 



-<5 3 (r) 



1 



3 WlqlTlQ 



((T q ■ <T Q ) 



(9) 



where Fq = ^ for quarks and F| = — A- for antiquarks. 
TTtq (toq) is the light (heavy) quark mass. The linear 
confinement potential is 



V^ nf = -4F c r F c n (a c „ n r+a 



cO \ 

aQ 1 T UqQ) 



There are similar expressions for and V^q . 

From Refs. [H, [Hij], one gets 



V»"( rij ) 



-C(g ch , m a , A)X 1 (m a ,A, nj), 



(10) 



Tfl 

C(g c h, m 7Ta , A) — — — X 2 (m na , A, r tj ) 
Y2m qi m qj 



x[<r(i) ■ o-(jWa(i)T a {j)], 



C(g chv ,m Pa ,A) 1 X 1 (m Pa ,A,r ij ) + ^ — 

I \)Tftq i Tflqj 



X 1 + 



fchv m qi + m qj ^ ^fchv ^ 2 m qi m qj 
9chv Mn g c hv 

xX2(m Pa , A, nj [r a (i)T a (j)}, 



(12) 



v uu( r ij) = C(9chv,m u ,,A)ixi(mu > A,r i j) + 

x (\ -\- ^ chv '^ m " a _|_ ( f chv ^ 2 m " 
V gchv Mn g C hv Mff 

xX 2 (m UJ ,A,n j )[cr(i) ■ a (J)] j, 



6m? 



yM 

uu 



(13) 
(14) 



C. Bound state problem 

According to the quark cluster model, the wave func- 
tion of the two mesons system in coordinate space reads 



* = ^ a ^a)^b{£,b)x{Rab)Z{R c 



(19) 



where £a = r 2 — >*i and ^3=^4 — are the internal 
coordinates of clusters A and B respectively, Rab is the 
relative coordinate between the two clusters, and R cm 
is the center of mass coordinate of the system. iPa(£a), 
and iPb{$,b) are the wave functions of A, B and Z(R cm ) 
represents the center of mass motion wave function of the 
system in coordinate space. All of them are treated as 
Gaussian functions: 



Mi a) 
Mis) 

Z (R cm ) 



3/4 



TOrOA 3 / 4 x m , ,*a 



(20) 



(11) where ttia = hib 



mqmQ is the reduced mass for the 

m q +niQ 

cluster A or B and Mab = M a +Mb = 2(m q +rriQ) is the 
total mass of the two clusters. The universal oscillator 
frequency uj is associated with the width parameter b u of 
the up quark through 



1 



(21) 



The unknown relative orbital wave function x(Rab) 
is expanded to partial waves 



x(Rab) = 



OO 

E 

N 



1 



Rab 

■(Rab) = 5>47ri? AS (^ 



X L (Rab)Y lm (Rab) 

/J,ABUJ\ 3 / i 



i=l 

XiL(VABOjRABSi), 



(22) 
(23) 
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where Si (i—l, 2,- • • , N) are the generator coordinates, 
fiAB = ^( m q + to q) is the reduced mass of the two clus- 
ters and i l (x) is the modified spherical Bessel function of 
L order. The coefficients Cj are to be obtained by solving 
the Schrodinger equation. 

The RGM equation for the bound state problem reads 

N 

Y}H^ - EN%] Cj =0 (i = 1, • • • , N) (24) 

3=1 

where 

{ % } = / YIm{ ^ ] { Nl } Y LM(S 3 )dS t dS 3 , 

{ S ■ } = / { ? } nSi) n drfci (25) 

with 

*(S<) = ^(CA)^(^)x(-RAB,S i )Z( J R cro p6) 
x(jRab , 5i ) = (^) 3/4 e -^x BW (H AB - Si f. (2?) 

Here 0a (4>b) denotes the total wave function of the clus- 
ter A (B), which includes the radial and spin parts. 

By solving Eq. the energy E and the correspond- 

ing relative motion wave function of the system (c$) are 
obtained. From the energy E, it is easy to derive the 
binding energy E = M D o + M D ,o — Mx D ■ If E is neg- 
ative, the system would be unbound. 



In the determination, we treat m u , m s , u> (or b u ), 
m Qi 9u and g s as inputs. For the up and strange quark 
masses, we use the values given in the previous work 
[H M, H HH rn u = 313 MeV and m s = 470 MeV. 
The width parameter b u — 0.5 fm in the chiral SU(3) 
quark model while b u = 0.45 fm in the extended chiral 
SU(3) quark model. These values have been fitted to 
reproduce the masses of the ground state baryons, the 
binding energy of the deuteron and the NN and YN 
scattering observables. To see the effects of this param- 
eter, we also use a larger value for the width parame- 
ter b u = 0.6 fm. To investigate the heavy quark mass 
dependence, we take several typical values m c = 1430 
MeV [13, m c = 1550 MeV [H| and m c = 1870 MeV 
[5j|. For the coupling constants, we can use (g u ,g s ) = 
(0.886,0.917), (0.880.755), (0.875,0.920), (0.237,0.451) 
or (0.363,0.500) (H,^. With these inputs, one gets sets 
of fitted values. Selected results are presented in Table [U 



m c (MeV) b 


- (fm) 


9u 


9c 


a c uc (MeV") 


oSi (MeV) 


1430 


0.45 


0.237 0.718 


45548 


-166.07 




0.5 


0.886 


0.774 


51320 


-143.36 




0.6 


0.886 


1.086 


54343 


-150.96 


1550 


0.45 


0.886 


0.642 


43129 


-150.08 




0.5 


0.886 


0.772 


47360 


-152.71 




0.6 


0.875 


1.097 


51600 


-161.86 


1870 


0.45 


0.363 


0.771 


27296 


-187.07 




0.5 


0.886 


0.858 


35445 


-168.59 




0.6 


0.886 


1.165 


42579 


-180.69 



TABLE I: Fitted parameters for the calculation in the hidden 
charm case. 



III. DETERMINING THE PARAMETERS 

There are numerous parameters in the Hamiltonian 
and the wave functions: g q , gq, a^g, a^q, WIq, m u , u, 
gch, gchv, fchv and A. The mass of the phenomenological 
a meson is also treated as an adjustable parameter. 

One should note, M D o, M D .o and Mx D are all calcu- 
lated with the Gaussian functions presented in the former 
section. The binding energy will be irrelevant to the in- 
ternal potentials of the color-singlet mesons because of 
the cancellation. That is, the form of the confinement 
and the values of g q , gQ, a c q Q and a^g will not give ef- 
fects to the numerical result of E . This feature can be 
understood with the effective potential between the clus- 
ters A and B in the generator coordinate method: 

V L 

V L (S i ,S j ) = -fi.-V D o-V D *o. (28) 

ij 

One can examine that the parts due to V OGE and V con f 
of Eq. © are exactly zero. Therefore, we may take 
any values, in principle, for these four parameters. In 
the following calculation, we deduce gQ, a c q Q and a c q q by 
fitting the masses of the ground state mesons D, D*, 
D s , D*, J/ip and rj c using a least square fit with the 
assumption a c cu — a c cs — a c cc . 



Actually, in the two meson molecule picture (see 
Fig.l), the meson exchanges play the dominant role in 
the energy of the system. The parameters of this part in- 
clude the quark-meson coupling constants and the meson 
masses. In the chiral quark model, the tt and a exchanges 
have the same coupling constant, named g c h, because of 
the chiral symmetry requirement. The coupling constant 
g c h is fixed through 

9ch _ 9 gjjNn ml, r^g) 
4tt 25 4tt m% 1 ' 

with gfj Nv / (Air) = 13.67 determined experimentally, 
from which one has g c h — 2.621. Thus, when the vector 
meson exchanges are not included,the mass of a and the 
cutoff mass A are the only adjustable parameters which 
can be fixed in the light quark systems. For the cou- 
pling constants of the vector meson exchange, one can 
use (g chv J chv ) =£3.0,0.0) 0, (2.09,5.26), (2.351,0.0), 
and (1.972,1.315) g|. 

To study the effects due to the uncertainty of the mass 
of the cr, we use m CT =595 MeV, 535 MeV and 547 MeV 
[54l [55|. For other mesons, we take the masses from the 
particle data book Q: m n o = 134.98 MeV, m p o = 775.8 
MeV and m w = 782.59 MeV. In the calculation, we adopt 
the cutoff A = 1000 MeV, 1100 MeV and 1500 MeV. 
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In this work, we also calculated the binding energy for 
the case of bottom analog where the flavor wave function 
is 



\Xb) = 



l r 

V2 



\B + B* 



\B* + B- 



(30) 



The procedure to determine the parameters is very simi- 
lar. Now the ground state mesons B, B* , B 3 , B*, T(1S) 
and rjh are used. As an input, we choose m^— 4720 MeV 
which is close to the value in Ref. [62j, m&=5100 MeV 
[H[ and m b =5259 MeV [5g[. By repeating the fitting 
procedure, one gets sets of parameters. We present the 
selected results in Table ITU 



m b (MeV) 6 


u (fm) 


9u 


9b 


a c ub (MeV^) 


aZ (MeV) 


4720 


0.45 


0.886 


0.897 


63452 


-149.44 




0.5 


0.875 


1.100 


68609 


-153.70 




0.6 


0.237 


1.567 


66134 


-207.04 


5100 


0.45 


0.875 


0.931 


46250 


-177.68 




0.5 


0.363 


1.208 


47881 


-200.74 




0.6 


0.875 


1.551 


61739 


-194.71 


5259 


0.45 


0.886 


0.943 


39073 


-183.09 




0.5 


0.363 


1.218 


41301 


-210.85 




0.6 


0.886 


1.565 


56580 


-204.13 



TABLE II: Fitted parameters for the calculation in the hidden 
bottom case. 
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FIG. 2: The effective potential V(S) for different meson ex- 
changes. The parameters used are b u — 0.5 fm, m u = 313 
MeV, m c = 1870 MeV, g ch = 2.621, g chv = 3.0, f chv = 0.0, 
m CT = 595 MeV, rn^ = 134.98 MeV, rn p = 775.8 MeV, 
= 782.59 MeV and A = 1100 MeV. Here the line for 
7r corresponds to V W (S), the line for a corresponds to V(S), 
and so on. 



IV. NUMERICAL RESULTS 

Before the numerical evaluation, we first take a look at 
the effective potential 

V(S) = V L =°(S, S) (31) 

where the generator coordinate S can qualitatively de- 
scribe the distance between the two clusters. These po- 
tentials rely on the meson exchange part in Eq. (jHJ. We 
illustrate the potentials corresponding to various consid- 
erations in Fig. [21 

From Fig [21 one sees that the interactions from tt, er, 
p and uj are all attractive [64]. The amplitudes for p 
and ui exchanges are comparable and their contributions 
should not be ignored arbitrarily. We will consider tt 
and a interactions for the moment and then include the 
vector meson contributions. 

Now we calculate the binding energy for the system 
D°D*° through solving Eq. (|2"4")) . Here we do not con- 
strain the sets of the parameters with the experimental 
data like the studies in Ref. [EU, [5||. After exploring 
all possible combinations of the parameters in the former 
section, we fail to get a bound state solution. Thus the 
D°D*° system is unbound when we consider only S-wave 
tt and a exchange interactions in this framework. This 
conclusion agrees with that of Ref. [4!| . 

Since the bottom quark is much heavier, the possibility 
of getting a bound state in the B meson systems is in- 
creased. Our former dynamical calculation is in favor of 



the existence of an S-wave Xb molecule. We also study 
this case in the present framework. When applying the 
evaluation to the bottom analog BB*, we get positive 
binding energies with the parameters in Section [Till The 
results are given in Table [TITl From that table, one finds 
a larger binding energy can be obtained with a larger m^, 
a smaller m a , a smaller b u or a bigger cutoff A. A deeper 
bound state should have a smaller root-mean-square ra- 
dius r rms , which is also illustrated in Table Hm 

To explore additional effects, we move on to include 
the vector meson exchanges. We use the parameters to 
reproduce experimental data for light quark systems [54| . 
The parameters and the results for D°D*° and B + B*~ 
are presented in Table IIVI and [V] respectively. For com- 
parison, the solutions without considering vector mesons 
are also given. Now a bound state seems to be possible 
in the D°D*° system. For its bottom analogy, the vector 
meson exchange interactions increase the binding energy 
about 10-20 MeV. 

Up till now, we considered only neutral components of 
the system. In Refs. (65| and [66[, the authors studied the 
case with symmetric wave function case (i.e. ao = a\ — 
-j= in Eq. [S]) and they found the coupling to charged 
components is important. We also present the numerical 
results for this case in Table IVII and I VIII which support 
the result that the channel coupling should be considered 
in studying X(3872). 
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m b (MeV) b u (fm) 


m CT (MeV) 


E (MeV) r rms (fm) 


4720 0.45 


595 


3.3/3.7/5.0 1.1/1.1/1.1 




547 


5.0/5.5/7.0 1.1/1.1/1.0 




535 


5.4/6.0/7.5 1.1/1.1/1.0 


0.5 


595 


2.0/2.3/3.1 1.3/1.3/1.2 




547 


3.4/3.8/4.7 1.2/1.2/1.2 




535 


3.8/4.2/5.2 1.2/1.2/1.2 


0.6 


595 


0.5/0.7/1.0 1.6/1.6/1.5 




547 


1.5/1.7/2.1 1.5/1.5/1.5 




535 


1.8/1.9/2.4 1.5/1.5/1.4 


5100 0.45 


595 


4.2/4.7/6.1 1.1/1.1/1.0 




547 


6.0/6.6/8.2 1.0/1.0/1.0 




535 


6.5/7.1/8.8 1.0/1.0/1.0 


0.5 


595 


2.8/3.1/4.0 1.2/1.2/1.2 




547 


4.2/4.6/5.7 1.2/1.2/1.1 




535 


4.6/5.1/6.2 1.2/1.1/1.1 


0.6 


595 


1.0/1.2/1.6 1.5/1.5/1.4 




547 


2.0/2.2/2.7 1.4/1.4/1.4 




535 


2.3/2.5/3.0 1.4/1.4/1.4 


5259 0.45 


595 


4.6/5.1/6.6 1.1/1.0/1.0 




547 


6.4/7.0/8.7 1.0/1.0/1.0 




535 


6.9/7.6/9.3 1.0/1.0/1.0 


0.5 


595 


3.1/3.4/4.4 1.2/1.2/1.1 




547 


4.6/5.0/6.1 1.1/1.1/1.1 




535 


5.0/5.4/6.5 1.1/1.1/1.1 


0.6 


595 


1.2/1.4/1.8 1.5/1.5/1.4 




547 


2.2/2.4/2.9 1.4/1.4/1.3 




535 


2.5/2.7/3.3 1.4/1.4/1.3 



TABLE III: Numerical results for the hidden bottom case 
when 7r and a exchange potentials are considered. The first 
(second,third) value for Eo and r rms corresponds to the cutoff 
A=1000 (1100,1500) MeV. 



X QM Ex. xQM 







fchv/ Qchv fchv J 


'g c hv = 2/3 


b u (fm) 


0.5 


0.45 


0.45 


m a (MeV) 


595 


535 


547 


Qchv 




2.351 


1.972 


E (MeV) 


X 


x/x/1.7 


X 


Trms (fm) 


X 


x/x/1.7 


x 


E (MeV) 


X 


x/0.5/3.0 


x 


r m s (fm) 


X 


x/1.8/1.6 


X 



TABLE IV: Numerical results for the hidden charm case when 
the vector meson exchange interactions are also included. 
Here m u = 313 MeV, g c h = 2.621 are used. The three values 
for Eq and r rm s correspond to m c = 1430 MeV, 1550 MeV 
and 1870 MeV in order. The first (second) Eo and r rms cor- 
respond to A=1100 (1500) MeV. x indicates the system is 
unbound. 



V. SUMMARY AND DISCUSSIONS 

In this work we have studied whether D°D*° (D°D*°) 
may form an S-wave molecule bound by the 7r, a, p and 
uj exchange interactions in a chiral quark model. These 
potentials are all attractive. By solving the RGM equa- 
tion, we failed to get a binding solution in this system 
if we consider only ir and a contributions. When the 
vector meson contributions are included, the existence of 



XQM Ex. xQM 

fchv I Qchv — fchv I Qchv — 2/3 

b u (fm) 05 045 045 

m CT (MeV) 595 535 547 

Qchv 2.351 1.972 

E (MeV) 2.3/3.1/3.4 17.2/18.7/19.3 11.1/12.3/12.8 

r rms (fm) 1.3/1.2/1.2 0.9/0.9/0.9 1.1/1.0/0.9 

E (MeV) 3.1/4.0/4.4 20.7/22.4/23.1 13.6/15.0/15.6 

r rmB (fm) 1.2/1.2/1.1 0.9/0.8/0.8 1.0/0.9/0.9 



TABLE V: Numerical results for the hidden bottom case when 
the vector meson exchange interactions are also included. 
Here m u = 313 MeV, g c h = 2.621 are used. The three values 
for Eo and r rms correspond to mi = 4720 MeV, 5100 MeV and 
5259 MeV in order. The first (second) Eo and r rms correspond 
to A=1100 (1500) MeV. 





xQM 


Ex. X QM 






fchv j Qchv 


fchvj Qchv — 2/3 


b u (fm) 


0.5 


0.45 


0.45 


m a (MeV) 


595 


535 


547 


Qchv 




2.351 


1.972 


Eo (MeV) 


X 


12.1/14.2/19.3 


4.5/6.0/9.7 


r rms (fm) 


X 


1.3/1.3/1.1 


1.6/1.5/1.3 


Eo (MeV) 


X 


16.3/18.6/24.5 


6.9/8.6/13.0 


r rms (fm) 


X 


1.2/1.2/1.0 


1.5/1.4/1.2 



TABLE VI: Numerical results for the hidden charm case with 
the symmetric wave function. Here m u — 313 MeV, g c h — 
2.621 are used. The three values for Eo and r rms correspond 
to m c = 1430 MeV, 1550 MeV and 1870 MeV in order. The 
first (second) set of Eo and r rms corresponds to A=1100 (1500) 
MeV. x indicates the system is unbound. 



D D molecule seems to be possible. The coupling to 
charged components is also important for a bound state. 

When moving on to the heavier B meson system, we 
obtain binding state solutions. Our calculation favors 
the existence of an S-wave BB* (BB*) molecular state, 
which agrees with the conclusion from Ref. [4!|. It will 
be very interesting to search for such a bound state in the 
radiative decay channel Xb — * B + B~-f and the strong 
decay channel Xb — > it + tt~T in the future. Finding 



XQM Ex. xQM 

fchvj ' Qchv = fchvj ' Qchv ~ 2/3 

bu (fm) 05 045 045 

m CT (MeV) 595 535 547 

gchv 2.351 1.972 

Eo (MeV) 12.4/13.8/14.4 47.6/50.0/50.9 32.0/34.0/34.8 

r rms (fm) 0.9/0.9/0.9 0.7/0.7/0.7 0.8/0.8/0.8 

Eq (MeV) 14.6/16.2/16.8 56.4/59.0/60.0 38.3/40.5/41.4 

r rma (fm) 0.9/0.9/0.8 0.7/0.6/0.6 0.7/0.7/0.7 



TABLE VII: Numerical results for the hidden bottom case 
with the symmetric wave function. Here m u = 313 MeV, 
g c h = 2.621 are used. The three values for Eo and r rma corre- 
spond to m b = 4720 MeV, 5100 MeV and 5259 MeV in order. 
The first (second) set of Eo and r rms corresponds to A=1100 
(1500) MeV. 
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it may be possible at the Tevatron or with the Large 
Hadron Collider beauty (LHCb) experiment [67l |. 

In the study of the deuteron, it was found that the ten- 
sor force which mixes the S-wave and D-wave interactions 
is crucial in binding the proton and the neutron. In an 
earlier calculation it was also concluded that the tensor 
potential is very important in the mesonic case (39j . In 
the present work, we did not consider effects from the D- 
wave. Further study using the current approach will be 
helpful to clarify whether this part can lead to a loosely 
bound D°D*° (D°D*°) state. 

From the numerical values, we observe that vector me- 
son contributions are important in binding two color- 
singlet mesons. However, the results rely on the vector 
coupling constants g c hv and f c hv Here we would like to 
mention that in our calculation the parameters of light 
quark part are taken from Ref. [54[, in which the calcu- 
lated NN scattering phase shifts and the binding energy 
of deuteron are consistent with the experimental data. 
But since the mechanism of the short range quark-quark 
interaction is still an open problem, whether OGE or 
vector meson exchange is dominate, or whether both of 
them are needed, one should be cautious when making 
conclusions from these results. 



In short summary, we have performed a dynamical cal- 
culation to investigate whether the D°D*° (D°D*°) may 
form a molecule by considering the n, cr, p and u> ex- 
change interactions. We could not find an S-wave molec- 
ular state in this system in the chiral quark model while 
its existence is not excluded in the extended chiral quark 
model. More details of the dynamics should be consid- 
ered in further study of the X(3872). If it is really not 
a molecule, the scheme of mixing a charmonium and a 
molecular state is probably a way to solve the puzzles of 
theX(3872). 
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